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Abstract 

Chemical turbulent mixing induced by rotation can affect the internal distribution of fi near the 
energy-generating core of main-sequence stars, having an effect on the evolutionary tracks 
similar to that of overshooting. However, this mixing also leads to a smoother chemical 
composition profile near the edge of the convective core, which is reflected in the behavior 
of the buoyancy frequency and, therefore, in the frequencies of gravity modes. We show that 
for rotational velocities typical of main-sequence B-type pulsating stars, the signature of a 
rotationally induced mixing significantly perturbs the spectrum of gravity modes and mixed 
modes, and can be distinguished from that of overshooting. The cases of high-order gravity 
modes in Slowly Pulsating B stars and of low-order g modes and mixed modes in j3 Cephei 
stars are discussed. 

Session: What can asteroseismology do to solve the problems 

Introduction 

Asteroseismology of main-sequence B-type stars is now providing us with constraints on the 
structure and on the internal rotation rate of massive stars (see e.g. the recent review by 
Aerts 2008). In particular, evidence for non-rigid rotation and for extra mixing near the core 
was found in several fi Cep pulsators (see e.g. Aerts et al. 2003, Mazumdaret al. 2006, 
Briquet et al. 2008, Dziembowski & Pamyatnykh 2008). In this context, an even deeper 
insight into the internal structure of B stars could be obtained if the seismic probe provided 
by the observed oscillation modes is fine enough to differentiate the near-core extra mixing 
due to overshoot from that induced by rotation. 

With the aim of studying the effects of different extra-mixing processes on the frequency 
spectra of SPB and fi Cep pulsators we focus, respectively, on models of 6 and 10 Mq. We 
consider the effect of overshooting and of turbulent mixing near the core that could be induced 
by the effects of rotation on the transport of chemicals. The models with overshooting were 
computed with CLES (Code Liegeois d'Evolution Stellaire, Scuflaire et al. 2008). In this code 
the thickness of the overshooting layer Aqv is parameterized in terms of the local pressure 
scale height H p : Aqv = «ov x (niin(r cc , H p (r cc )), where r cc is the radius of the convective 
core given by the Schwarzschild criterion and c*ov is a free parameter. In the overshooting 
region mixing is assumed to be instantaneous and the temperature stratification radiative. 

As a first approach, the turbulent mixing has been modelled in CLES by a diffusion pro- 
cess with a parametric turbulent diffusion coefficient Dx that is uniform inside the star and 



2 Discriminating between overshooting and rotational mixing in massive stars: any help from asteroseismology? 




Figure 1: Left panel: Hydrogen mass fraction profile in the central regions of 10 Mq models with Xq ~ 
0.3. The different lines correspond to models computed with CLES, with overshooting (dashed line), and 
turbulent diffusion coefficient Dt (solid line), and with the Geneva code with an initial rotational velocity 
of 50 km s -1 (dotted line). Right panel: Brunt- Vaisala frequency profile in the central regions of the 
10 Mq CLES models shown in the left panel. 



independent of age. We computed also models using the Geneva code that includes the treat- 
ment of rotation and related transport processes as described in Eggenberger et al. 2008. The 
comparison between the two series of models shows that the chemical composition profiles in 
the central regions provided by a uniform and time independent diffusion coefficient represent 
a first approximation, at least for massive stars, of the effect that would be produced by 
such a rotationally induced chemical transport (see Fig. 1). The values of the parameter 
Z>T were chosen in order to be close to the value of the chemical diffusion coefficient near 
the core provided by the Geneva models (see Montalban et al. 2008 for more details on this 
parametrization). 

SPBs and /3 Cep pulsators are considered as slow or moderate rotators. SPBs show a 
typical rotational velocity of 25 km s _1 (Briquet et al. 2007), whereas the range of projected 
rotational velocity in /3 Cep stars extends from to 300 km s _1 with an average of 100 km s _1 
(Stankov & Handler 2005). The Geneva code calculations for 10 Mq models provide values 
of the chemical diffusion coefficient near the convective core with Xq = 0.3, of the order of 
5 X 10 4 cm 2 s — 1 for a rotational velocity on the zero-age main sequence (hereafter named 
initial rotational velocity, V;) of 20 km s _1 , 7 X 10 4 cm 2 s _1 for Vi=50 km s _1 , and 1.6 X 
10 5 cm 2 s — 1 for Vi=100 km s _1 . On the other hand, the effect of an initial rotational velocity 
of 25km s _1 on the central hydrogen distribution of a 6 Mq model is well mimicked by a 
Dt ~ 5000 cm 2 s _1 . The results presented in this paper concern mainly the parametric 
models with D T = 5 X 10 3 cm 2 s- 1 for SPB models and D T = 7 X 10 4 cm 2 s- 1 for fi Cep 
ones, and they will be compared with overshooting models closely located in the HR diagram, 
that means arjv=0 for the SPB case and 0.1 for the fi Cep one. 

As shown in Fig. 1, models located at the same position in the HR diagram, but computed 
with overshooting or with turbulent mixing, show a significantly different chemical composi- 
tion profile near the core (see Fig. 1 left panel) and consequently a different behaviour of the 
sharp feature in the Brunt- Vaisala frequency (N) located in the ^-gradient region (Fig. 1 
right panel). 
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Figure 2: Behaviour of the Brunt- Vaisala frequency (left panels), of the hydrogen abundance profile (middle 
panels) and of the £ — 1 g-mode period spacing in models of 6 Mq (right panels). We consider several 
models along the main sequence with decreasing central hydrogen abundance {X c — 0.7, 0.5, 0.3 and 0.1). 
In the right panels the uppermost model has X c — 0.7 and the lowermost X c — 0.1. Upper panels refer 
to models computed with no extra mixing, whereas lower panels in models of 6.0 MQComputed with a 
turbulent diffusion coefficient Dt — 5000 cm 2 s — 1 . 



High-order g modes in Slowly Pulsating B stars 

As presented in Miglio et al. (2008a) (hereafter Paper I), the periods of high-order gravity 
modes in main-sequence stars can be related, by means of analytical expressions, to the 
detailed characteristics of the /x-gradient region that develops near the energy-generating 
core, and thus to the mixing processes that affect the behaviour of fi in the central regions. 
The period spacing (AP) of high-order g modes can be described as a superposition of a 
constant term predicted by the first order approximation by Tassoul (1980) and periodic 
components directly related to the location and sharpness of V M . 

We recall that in the asymptotic approximation presented in Tassoul (1980) the periods 
of high-order gravity modes are given by P^i = 7r 2 IIo L~ 1 (2k + cte), where L = + 1)] 1 / 2 
(with I the mode degree), k the order of the mode and II^ 1 = J 1 -p-dx' , where x' is the 
normalized radius and xq corresponds to the boundary of the convective core. 

In main-sequence stars V/i depends on the size of the convective core and on extra- 
mixing processes (e.g. overshooting, diffusion, turbulent mixing) that may alter the chemical 
composition profile in the central region of the star. We here briefly describe the effect 
of rotationally induced mixing by considering models of a 6 Mq star computed including 
turbulent diffusion near the core (described by a turbulent-diffusion coefficient Dx = 5 X 10 3 
cm 2 s~ 1 ), which mimics the effects of an initial rotational velocity of 25km s — 1 . In Fig. 2 
we compare the period spacing of such models (lower-right panel) with models that share the 
same location in the HR diagram, but are computed without any extra-mixing (upper-right 
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Figure 3: Frequencies of pulsation modes with angular degree t — — 3 as a function of log T c ff for 
main-sequence models of a 10 Mq star. Gray dots correspond to the frequencies of models computed 
with an overshooting parameter aov = 0.1, whereas black dots are the frequencies of models computed 
with a turbulent diffusion coefficient D T = 7 X 10 4 cm s~ 2 . Excited modes are represented by thicker 
symbols. The vertical lines indicate the effective temperature of models with a hydrogen mass fraction at 
the center of the order of 0.5, 0.3 and 0.1 (left to right). 



panel). We notice that turbulent mixing has a substantial effect on the period spacing: the 
amplitude of the periodic components in AP becomes a decreasing function of the radial 
order k. This behaviour can be easily explained by the analytical approximation in Paper I, 
provided that the smoother fi profile (and thus N, see Fig. 2) is taken into account in the 
analysis. 



Low-order g modes and mixed modes in f3 Cephei stars 

As recalled in the previous section, changes of fi induced by the turbulent mixing acting 
near the core lead to significant variations in the properties of the spectrum of high-order g 
modes. Although the approximation presented in Paper I and used to relate the sharpness of 
the "bump" in the Brunt- Vaisala frequency to the properties of g modes is no longer valid 
for low-order modes, the analytical description of gravity modes mentioned there is still able 
to qualitatively relate some properties of low-order g-mode spectra to the characteristics of 
the /^-gradient region near the core (see Miglio et al. 2008b). 
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Figure 4: Frequency difference (in c/d) between 
modes with consecutive radial order (k) and de- 
gree I = 2 for the 10 Mq cles models in Fig. 1. 
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During the main sequence, the combined action of nuclear reactions and convective mixing 
leads to a chemical composition gradient at the boundary of the convective core, to a decrease 
of Iln, and thus to an increase of the frequencies of gravity modes. The latter interact with 
pressure modes of similar frequency and affect the properties of non-radial oscillations by the 
so-called avoided-crossing phenomenon. The modes undergoing an avoided crossing (mixed 
modes) are therefore sensitive probes of the core structure of the star. 

As presented in Montalban et al. (2008), we describe the effects of extra-mixing on the 
oscillation spectrum of a typical f3 Cep star by comparing the properties of low-order g and p 
modes in models computed with overshooting and with chemical turbulent diffusion. Since 
the parameters aov ar| d £*T were chosen to lead to similar evolutionary tracks on the HR 
diagram, such comparison allows us to remove the differences in the frequencies due to a 
different stellar radius. In Fig. 3 we plot the oscillation frequencies for 10 Mq models along 
the main-sequence phase. As expected, the differences between frequencies of overshooting 
models and turbulent mixing ones is very small for pressure modes, while significant differences 
appear for gravity and mixed modes. These differences increase with the radial order of g 
modes and, for a given frequency, with the angular degree of the modes. In Fig. 3 we have 
also marked with thicker symbols the frequencies of modes that the non-adiabatic oscillation 
code MAD (Dupret et al. 2003) predicts to be excited. Here we show only the case with 
£>t = 7 X 10 4 cm 2 s _1 but computations with lower or higher efficiency of the turbulent 
mixing show that the differences between overshooting and turbulent-mixing models increase 
with the value of D-p. 

What is most relevant from an asteroseismic point of view, is the change in the distance 
between consecutive frequencies of g-modes or modes of mixed p-g character (Aj/). In 
Fig. 4 we plot these differences for the 1 = 2 modes of the 10 Mq models with an effective 
temperature T c ff ~ 22550 K (Xq ~ 0.3). The dots indicate differences computed between 
pairs of excited modes. Therefore the difference of Au that we can expect between models 
with sharp V M (overshooting models for instance) and models with a chemical composition 
gradient smoothed by the effect of, for instance, a slow rotation (V ro t ~ 50 km/s) is of 
the order of 0.4 c/d (~ 5fiHz), much larger than the precision of present and forthcoming 
observations. 
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Conclusions 

We presented how the frequencies of gravity modes and mixed modes in ft Cep and SPB stars 
depend on the detailed characteristics of the ^t-gradient region that develops near the energy 
generating core, and thus on the mixing processes that can affect the behaviour of /i in the 
central regions. In particular we have shown that for rotational velocities typical of main- 
sequence B-type pulsating stars, the signature of a rotationally induced mixing significantly 
perturbs the spectrum of gravity modes and mixed modes, and can be distinguished from 
that of overshooting. 

Such a sensitivity of g-mode oscillation frequencies to near-core mixing can provide an 
additional constraint to the modelling of massive stars with rotation, especially when coupled 
with seismic inferences on the internal rotational profile (see e.g. the review by Goupil &i 
Talon 2009), and with constraints on chemical enrichments in the photosphere (e.g. Maeder 
et al. 2008 and Morel 2008). Further investigations are however needed to assess under which 
observational conditions such information can be recovered from the oscillation frequencies, 
given realistic uncertainties on stellar global parameters, errors on oscillation frequencies and, 
in the case of SPBs, given the severe influence of rotation on the spectrum of g modes. 
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